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We examine the quantum Hall effect in bilayer graphene grown on Cu substrates by chemical 
vapor deposition. Spatially resolved Raman spectroscopy suggest a mixture of Bernal (A-B) stacked 
and rotationally faulted (twisted) domains. Magnetotransport measurements performed on bilayer 
domains with a wide 2D band reveal quantum Hall states (QHSs) at filling factors v — 4, 8, 12 
consistent with a Bernal stacked bilayer, while magnetotransport measurements in bilayer domains 
defined by a narrow 2D band show a superposition of QHSs of two independent monolayers. The 
analysis of the Shubnikov-de Haas oscillations measured in twisted graphene bilayers provides the 
carrier density in each layer as a function of the gate bias and the inter-layer capacitance. 

PACS numbers: 73.22.Pr, 73.43.-f, 68.35.bp 



Bilayer graphene with Bernal-type stacking is an in- 
teresting electron system possessing high carrier mobil- 
ity and tunable energy band-gap, 1 4 which makes it a 
potential material for electronic devices. Recent studies 
showed large area growth of bilayer graphene on SiC and 
metal substrates by chemical vapor deposition (CVD). 
The growth of Bernal stacked bilayer graphene on SiC 
substrates has been confirmed by electron microscopy^ 
and magneto-transport measurements. 6 While several re- 
cent studies suggest the growth of Bernal stacked bilayer 
on metal substrates based on Raman spectroscopy, 7 10 
evidence of stacking from electron transport data in 
grown bilayer graphene has been scant. It is therefore 
interesting to probe the electronic properties of CVD- 
grown graphene bilayers, which in turn can shed light on 
the growth mechanism and help assess its potential for 
device applications. Here we provide a systematic inves- 
tigation of the quantum Hall effect in bilayer graphene 
grown on Cu substrates by chemical vapor deposition. 
Our data show that such bilayers consist of a mixture of 
domains which are either Bernal stacked or are rotation- 
ally faulted ('twisted bilayer'). 

The graphene samples studied here are grown on a 25 
/im-thick Cu foil at a temperature of 1035 °C by CVD, 
using a mixture of methane and hydrogen at the partial 
pressures of 0.02 mbar and 0.03 mbar, respectively. Af- 
ter the growth, the graphene film on one side of the Cu 
foil is coated with PMMA and placed in an aqueous so- 
lution of ammonia persulfate ((NH4)2S2 08) to dissolve 
the Cu on the unprotected side. The PMMA film that 
carries the graphene flake is rinsed several times with de- 
ionized water to minimize the chemical contamination, 
and then transferred onto a silicon substrate covered with 
285 nm-thick thermally grown Si02- After the transfer, 
the sample is allowed to dry and the PMMA is dissolved 
in acetone. 

Figure 1(a) shows the optical micrograph of the 
graphene film transferred on the Si02 substrate. The 
optical inspection indicates the presence of monolayer, 



bilayer, and trilayer regions. To probe the number of 
graphene layers, and obtain an initial assessment of the 
layer stacking, the sample is characterized by Raman 
spectroscopy acquired using a 488 nm excitation wave- 
length, 300 nm spot size, and a power lower than 0.1 mW. 
Figure 1(b) presents a representative mapping of the Ra- 
man 2D band full width at half maximum (FWHM) ac- 
quired over a 30x30 jam 2 area. These data reveal the 
presence of distinct domains on the bilayer area with ei- 
ther a narrow 2D band, with FWHM values between 27 
and 33 cm -1 , or a wide 2D band, with FWHM values 
between 45 and 54 cm -1 . By comparison the Raman 
2D FWHM measured in monolayer graphene is 28 — 30 
cm -1 while in Bernal (A-B) stacked bilayer graphene the 
2D band is wider ilr— Figure 1(b) data therefore sug- 
gest that the bilayer domains with narrow 2D band con- 
sist of two graphene monolayers which are rotationally 
faulted (twisted bilayer), while the domains character- 
ized by wider 2D band consist of two Bernal stacked 
monolayers. We note the two types of bilayer domains of 
Fig. 1(b) show no obvious differences in optical contrast. 
Figure 1(c) shows samples of Raman spectra acquired 
on the same sample of Fig. 1(b), at different positions 
on the monolayer, the twisted bilayer, and the Bernal 
stacked bilayer regions, as indicated. The 2D FWHM 
of these Raman spectra are 28 cm -1 , 30 cm -1 , and 50 
cm -1 , respectively. The 2D band intensity (ho) is larger 
than the G band intensity (Iq) on the monolayer and 
bilayer domains with narrow 2D band, an observation 
which agrees with Raman spectroscopy results for exfoli- 
ated graphene j n i 12 In contrast, the bilayer domain with 
a wide 2D band shows an I^dI^g ratio lower than 1. The 
D band, located at a Raman shift of 1350 cm -1 is either 
absent or very weak, indicating that the defect density of 
the CVD-grown and transferred graphene is low. 

Figure 1(d) shows a histogram of the 2D band FWHM 
values acquired over a 15x20 jim 2 bilayer graphene grain 
characterized by a wide 2D band. The data points range 
between 45 cm -1 and 65 cm -1 , with a maximum at 53 
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FIG. 1: (Color online) (a) Optical micrograph of a multi-layer graphene sample after transfer onto a 285-nm thick SiCb/Si 
substrate. The dashed lines delineate monolayer, bilayer, and trilayer regions, (b) 2D FWHM map of graphene over a 30x30 
fim 2 area showing the bilayer is a mixture of domains with either wide 2D band (45 — 54 cm -1 ) or narrow 2D band (27 — 33 
cm -1 ). The dotted line marks a Hall bar subsequently realized to probe electron transport in individual bilayer regions, (c) 
Raman spectra acquired at three different positions, as marked in panel (b) show the G (~ 1580 cm -1 ) and 2D (~ 2700 cm -1 ) 
bands, (d) Histogram of the 2D band FWHM on a bilayer domain with wide 2D band. The average 2D FWHM is 53 ± 2 
cm -1 , (e) Example of a 2D band spectrum (black line) acquired on a Bernal stacked bilayer domain. A fit (red) using four 
Lorentzian functions (green) provide a very good match to the experimental data, (f) Optical micrograph of a back-gated Hall 
bar fabricated on bilayer graphene. 



cm -1 . Figure 1(e) presents a typical spectrum of the 
2D band selected from the bilayer graphene region with 
a wide 2D band. The 2D band spectrum of Fig. 1(e) 
could not be fitted with a single Lorentzian function, but 
an excellent fit is obtained using four Lorentzian func- 
tions. The combined data of Fig. l(b-e) therefore sug- 
gest that bilayer domains with narrow 2D band consist of 
twisted graphene monolayers, while the bilayer domains 
with wide 2D band are two Bernal stacked monolayers. 
We next focus on the electron transport and quantum 
Hall effect in these two types of bilayer s. 

After the graphene is characterized by Raman spec- 
troscopy, we fabricate back-gated Hall bar devices on se- 
lected bilayer domains with a uniform 2D peak FWHM, 
which is either narrow (27 — 33 cm -1 ) or wide (45 — 65 
cm -1 ). The active region of the Hall bar is defined by 
electron-beam (e-beam) lithography and isolated from 
the rest of the film using oxygen plasma etching. Metal 
(Ni) contacts are defined by a second e-beam lithogra- 
phy, metal deposition, and lift-off. An example of a 
Hall bar device made on a bilayer graphene is shown in 
Fig. 1(f). The carrier mobility (/i) of each sample is 
determined from the four-point conductivity (a) depen- 
dence on back-gate bias (Vbg), M = 1/Cbg x da/dYsG] 



Cbg is the back-gate capacitance per unit area. For the 
samples examined in this study Cbg = 14.4 nF/cm 2 , a 
value which is measured on 100 x 100 fim 2 metal pads de- 
posited in close proximity of the Hall bar devices, and 
also confirmed by Hall measurements. The extracted 
mobility of the Bernal stacked bilayer graphene devices 
are in the 700 — 1, 800 cm 2 /Vs range, and twisted bilay- 
ers exhibit mobility values in the 3, 400 — 3, 700 cm 2 /Vs 
range. The higher mobility in twisted bilayer s by com- 
parison to Bernal stacked bilayers can be explained by 
differences in their band- structure, which forbid elec- 
tron back-scattering in monolayer graphene, and hence 
in twisted bilayer, but allow back-scattering in Bernal 
stacked bilyersJ^ 

To establish the layer stacking of the CVD grown 
graphene bilayers and explore their electronic properties, 
in the following we discuss quantum Hall effect measure- 
ments in this system. Magnetotransport measurements 
were carried out in perpendicular magnetic fields (B) up 
to 31 T, using a pumped 3 He refrigerator with a base 
temperature T = 0.3 K, and small signal, low frequency 
lock-in techniques. Figure 2(a) shows the longitudinal 
(pxx) and Hall resistivity (p X y) as a function of Vbg mea- 
sured at a perpendicular magnetic field B = 25 T and 
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FIG. 2: (Color online) (a) p xx and vs. Vbg, measured at 
B = 25 T, and T = 0.3 K. The data shows QHSs, marked 
by vanishing p xx and quantized p xy at filling factors v — ±4 
and v — —8. (b) p xx vs. v measured at T = 0.3 K, and at 
different B-field values, illustrating the emergence of QHSs 
at integer filling factors multiple of four with increasing the 
B-field. 



FIG. 3: (Color online) ( a ) p xx vs. Vbg measured at different 
temperatures, and at B = 25 T. Inset: p xx vs. T _1 at v — — 4 
on a log-lin scale, measured at B — 15, 20, 25, 30 T. The p xx 
temperature dependence is determined by the QHS energy 
gap (A), (b) A vs. B, for v = -4 and v = -8 QHSs. 



T = 0.3 K in a graphene bilayer that displays a Raman 
signature consistent with Bernal stacking, i.e. wide 2D 
band. The data shows clear quantum Hall states (QHSs), 
marked by vanishing p xx and quantized p xy at filling fac- 
tors v = ±4 and v = — 8. The filling factors are deter- 
mined from the p xy plateau values, which are equal to 
h/ve 2 ; h is Planck's constant and e the electron charge. 
Alternatively, the filling factor can be calculated using 
v = nh/eB, where n is the total carrier density calcu- 
lated by n = Cbg(Vbg — Vb)/e; Vd is the gate bias at 
the charge neutrality (Dirac) point. 

Figure 2(b) shows the p xx vs. v measured in the same 
sample at different B values, and at T = 0.3 K. The data 
show the emergence of QHSs at integer filling factors that 
are multiples of four, i.e. v = ±4, —8, —12, thanks to the 
fourfold degeneracy of each Landau level (LL) associated 
with the spin and valley degrees of freedom^ The QHSs 
filling factors of Fig. 2 are consistent with the expected 
values in bilayer graphene ^ 17 ! 18 in effect fingerprinting a 
Bernal stacked bilayer. 

Figure 3(a) shows p xx vs. Vbg measured in the same 
sample discussed in Fig. 2 at B = 25 T, and at different 
temperatures. Although the v = ±4, —8 QHSs weaken 



with increasing T, these QHSs remain clearly visible at 
the highest temperature, T = 70 K. The inset of Fig. 3(a) 
shows the Arrhenius plot of p xx measured at v — — 4, and 
at B = 15,20,25,30 T. These data follow a thermally 
activated behavior, p xx oc e - A /( 2k BT) ^ w \ iere /\ [ s th e 
energy gap and ks is Boltzmann's constant. Figure 3(b) 
shows the extracted v = —4, —8 QHSs energy gaps vs. B. 
The data follow a linear dependence of A as a function 
of 5, with the A values approaching at B c± 13 T. The 
QHS energy gaps of Fig. 3(b) are considerably smaller 
than theoretical expected values. 19 For example, the the- 
oretically expected energy gap of v = — 4 at B = 30 T is 
108 meV, a value roughly eight times larger than the ex- 
perimental value. The QHS energy gaps probed in CVD- 
grown bilayer graphene are also smaller, albeit compara- 
ble than the values typically measured in exfoliated bi- 
layer graphene on Si02 substrates 3*21 The smaller QHS 
energy gap values can be explained by LL broadening due 
to disorder, most likely charged impurities in the prox- 
imity of graphene, and defects . 22 i 23 

We now turn to the magneto-transport properties of 
the twisted bilayer graphene samples, fabricated on bi- 
layer graphene domains with a narrow Raman 2D band. 
Figure 4(a) shows an example of p xx and p xy vs. B data, 
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FIG. 4: (Color online) (a) p xx and p xy vs. B measured at n = 9.7 x 10 12 cm -2 and at T = 0.3 K. The SdH oscillations stem 
from a QHSs superposition of the two decoupled graphene monolayers, (b) Fourier transform of p xx vs. B~ x data. The two 
peaks represent the layer densities, up to a factor Ae/h. (c) Top layer, bottom layer, and total carrier density of the twisted 
bilayer graphene vs. Vbg- The symbols are experimental data, and the solid lines are calculations. 



measured in a twisted bilayer device at Vbg — — 40 V, 
at a carrier concentration n = 9.7 x 10 12 cm -2 , and at 
T = 0.3 K. These data possess several noteworthy fea- 
tures. First, the p xx vs. B data display Shubnikov-de 
Haas (SdH) oscillations present down to magnetic fields 
as low as B c± 3 T, which contrast Figs. 2 and 3 data, 
where QHSs are not visible at 5-fields lower than 10 T. 
This contrast can be explained by the larger monolayer 
graphene LL energies by comparison to bilayer graphene. 
Moreover, the p xx vs. B data do not follow a QHS se- 
quence which can be readily attributed to either mono- 
layer (y = ±2,6,10...) or bilayer [y = ±4,8,12...), and 
instead shows a beating pattern similar to the QHSs su- 
perposition of a multisubband system. This observation 
is consistent with parallel electron transport in two inde- 
pendent graphene monolayers. 

To determine the subband (layer) density in a twisted 
bilayer, we examined the Fourier transform (FT) of p xx 
vs. B~ x data, calculated by first re-plotting the p xx vs. 
B data, subtracting a linear fit background to center the 
p xx vs. B~ x data around zero, and then applying a fast 
Fourier transform (FT) algorithm. Figure 4(b) shows the 
FT amplitude vs. B corresponding to Fig. 4(a) data. 
These data show two prominent peaks, which yield the 
two layer densities, up to a factor 4(e/ft) = 9.67 x 10 10 
cm _2 -T _1 . Figure 4(c) summarizes the layer and total 
densities as a function of Vbg- We attribute the higher 
(lower) density to the bottom (top) layer, as it lies closer 
(farther) with respect to the back-gate. Both layer den- 
sities go to zero at Vd = 68 V. 

To understand the layer density dependence on Vbg 
in twisted bilayer graphene, we employ a model used to 
calculate the layer densities in independently contacted 
graphene double layers separated by a dielectric. 24 The 
applied Vbg is distributed partly across the SiC>2 dielec- 
tric and partly on the Fermi energy of the bottom layer: 

e(V BG ~ V D ) = e\n T ± n B ))/C B G ± E F (n B ) (1) 

Here E F (n) = sgn(n)vF\^m is the Fermi energy relative 




FIG. 5: (Color online) Twisted bilayer p xx contour plot as a 
function of Vbg and B. The black (red) lines are the calcu- 
lated position of v — ±0,4,8,12... LLs of the bottom (top) 
layer. The parentheses show the corresponding filling factors 
of the bottom and top layers, respectively. 



to the charge neutrality point in monolayer graphene at a 
carrier density n; sgn represents the sign function. Sim- 
ilarly, the bottom layer Fermi energy is the sum of the 
electrostatic potential difference between the layer and 
the Fermi energy of the top layer: 

E F (n B ) = e 2 n T /C int ± E F (n T ) (2) 

where Cm is the interlayer capacitance. Using Eqs. (1, 
2) and Cm as a fitting parameter, we calculate n# and 
tit as a function of Vbg- An excellent fit to the exper- 
imental data is obtained for Cm = 6.9 /iF/cm 2 [solid 
lines in Fig. 4(c)]. Remarkably, this value is in good 
agreement with the inter-layer capacitance expected for 
a Bernal stacked bilayer, 25 suggesting that the separation 
of the two layers in twisted bilayer graphene is close to 
that of a Bernal stacked bilayer. A similar Cm value has 
been reported in a sample where two graphene monolay- 
ers were sequentially transferred onto each other ^ 
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Figure 5 shows the p xx contour plot as a function of 
Vbg an d B probed in the twisted bilayer sample of Fig. 4. 
The charge neutrality (Dirac) point is reached at back- 
gate bias Vd — 68 V. The data show a QHS pattern 
which stems from the QHSs superposition of the two 
decoupled monolayers. To quantitatively map the po- 
sition of the observed QHSs, we use Eqs. (1) and (2) to 
calculate the layers densities as a function of Vbg and 
B, with the only difference that the Fermi energy de- 
pends on both density and magnetic field as Ep = E^, 
where En = sgn(N)vF ^/2e%B\N\ is the energy of the 
N th LL in monolayer graphene, and TV = Int[nh/4eB] 
is the LL index; Int is the nearest integer function. Us- 
ing Ci n t =6.9 /iF/cm 2 , the value extracted from Fig. 4 
data analysis, at fixed B and Vbg values we calculate 
riB and nx, which are then converted into layer filling 
factors vb,t — riB^rh/eB. The black (red) lines in Fig. 
5 represent the calculated position of half-filled LLs, i.e. 
vb,t = ±0,4,8,12... for the bottom (top) layer. The 



p xx maxima are in excellent agreement with the calcula- 
tions, quantitatively confirming that the QHS sequence 
of twisted bilayer graphene is a superposition of the QHSs 
of the two graphene monolayers. 

In summary, using a combination of Raman spec- 
troscopy and magnetotransport measurements we estab- 
lished that CVD-grown bilayer graphene on Cu consists 
of a mixture of Bernal stacked and twisted monolayer 
domains. The Bernal stacked domains show QHSs at 
filling factors v = 4, 8, 12 in agreement with data in ex- 
foliated bilayer graphene. The twisted bilayer graphene 
domains display a superposition of the individual QHSs 
of two grapehene monolayers, which allows us to extract 
the layer densities and inter-layer capacitance. 
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